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О сходимости экстремальных методов 
в бесконечномерных задачах оптимизации 
с параболическим уравнением 


В статье рассматривается сходимость экстремальных методов в задаче оптимизации системы с квазилинейным 
параболическим уравнением. На примере бесконечномерной минимизации квадратичного целевого 
функционала показано, что традиционные методы, включающие методы наискорейшего спуска, 
МОНОТОННОГО убывания и сопряженных градиентов, сходятся не равномерно. Равномерной СХОДИМОСТИ 
позволяет добиться прямой экстремальный метод с регулируемым направлением спуска. 


Рассмотрим проблемы сходимости классических экстремальных методов на 
примере задачи оптимального управления тепловым потоком вдоль вертикального 
непрерывного цилиндрического слитка. Это одна их актуальных задач металлургии. 
Установившиеся тепловые процессы в таких задачах можно описать следующим 
квазилинейным параболическим т 


тт ив) =0, (7,2) ЕО =[0,^]х[ 0,7] (1) 
с граничными условиями: 
ОТ ОТ ОТ 
РЯ — 0, ети ПЕ и(2), = УТ Г), Т 0<„<в = То ; (2) 
ог 0, "|; "|; =—0 


Здесь 5 = {2 А Ей } — область определения управляемого теплового 
потока и(2), влияющего на качество слитка, 5, = Ах[ 0,2, |-— область интенсивного 
заданного (не управляемого) охлаждения слитка, Т — температура в слитке, А-— 
радиус слитка. __ остальные параметры уравнения (1) и граничных условий (2) — 
&(Т), В(Г), у, — считаем известными и корректно заданными [1]. 

Тестовую ен оптимального управления сформулируем как задачу синтеза 
заданной температуры Т,(7,2) в части слитка [0, К]х[2,,0]. Необходимо в граничном 
условии (2) найти функцию и(2), доставляющую минимум целевому функционалу 


КЕ 
1= | [(г-т.) анё. (3) 
02, 
Тестовые расчеты организовывались следующим образом. Задавалось гранич- 
ное управление и(2), которое считалось оптимальным их(2), и по уравнению (1) с 
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граничными условиями (2) рассчитывалась соответствующее Т,(г,2). Далее предпола- 
галось, что и»(2) неизвестно и для его определения необходимо решить бесконечно- 
мерную экстремальную задачу 
Л(и(2)) —> шщ. 

Для минимизации на открытых множествах управлений, обычно, используют 

итерационные градиентные методы [2-4]: 
и (2) =и" (2) -№“УЛи“:2) наб, 5“ >00, К=0,1,.... (4) 

В случае выбора параметра шага Ь^ из условия наилучшего продвижения в за- 

данном направлении минимизации — У.Л (и ;2), типа 


Ь" = аге той (и" — БУ. (и";2)), 


мы имеем метод наискорейшего спуска. 

В результате применения методики, изложенной в работе [5], применительно к 
целевому функционалу (3) был получен градиент У./ ‚ который определяется из вы- 
ражения: 

УЛ (и"; 2) =УЛ" =-р, (5) 
где сопряженная переменная й(’,2) удовлетворяет следующей сопряженной задаче: 


19 ой ой 
тата а 


(6) 
ОВ(Г) ОТ ой де(Т) Т тт.) =0, ео 
ОТ ди дк ОТ 02 
Граничные условия имеют вид 
И = и — и” =0::й|. = 0% (7) 
Оу |"=0 в 5 Я |; т” 


0<2<й 

Метод наискорейшего спуска (4) для минимизации функционала /(и) реализо- 
вывался приблизительно с условием: 

если Л“ <.Л\, тогда БН =ЬЬ, В >1; 

тогда повторяется предыдущая 8 
если Л" > Л*, ыы ь (8) 
итерация, при Б’ =Ь.Б`; О<Ь, <1. 
Принималось ВБ, = 1,2, Б, = 0,5, Ь° задавалось из условия не более 20 % первого изме- 
нения функции и(2) на первой итерации [5], [6]: 
0,2м° (2) 
У. (и°;2) 


0 
На рис. | показана сходимость функции и(2) от начального приближения и (2) = 


В, = шт 


2 


. (9) 


= 200 кДж/(м°с) к точному решению и, на 5, 15 и 100 итерациях. Видно, что градиент- 
ный метод наискорейшего спуска (его приблизительная реализация) (4), (8) не сходится к 
точному решению и, за значительное число итераций. От итерации к итерации схо- 
димость метода существенно замедляется. В табл. 1 приведены некоторые значения 


целевого функционала (3) и невязки решений Ли” = м ео 
221» 


Вряд ли следует ожидать существенного улучшения решения при дальнейшем 
увеличении числа итераций, к тому же разумным числом следует считать число ите- 
раций не более 100. После 100 итераций максимальное расхождение найденной и 


= 42 К. 


требуемой функции составляет тах/Т -— Т. 
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Часто в градиентных алгоритмах минимизации используют не метод наиско- 
рейшего спуска, а так называемый метод монотонного убывания [2-4], в котором 


принимают параметр 5“ = ° постоянным на всех итерациях. Как показывают расче- 
ты (табл. 1), такой метод дает существенно худшие результаты. Действительно, после 
100 итераций пах|Т —Т. 


= 63 К, что заметно хуже предыдущего результата. 


Не удалось, по существу, улучшить решение и применением метода сопряжен- 
ных градиентов 


и" (2)=и"(2)-В' р’ (2), (0 


где р*=УЛ МЛ ЛУ ре, р°=УЛ. 


м, «Дж с) 


щ =" 
125 чай 
мыть 


0 5 10 2,м 


Рисунок 1 — Сходимость функции и(2) к точному решению и. (2) 
при использовании метода (4), (8) 


Решение данным методом представлено в табл. | и на рис. 2. Здесь Ь, = 0,4. Кроме 
того, нижняя строка условия (8), для очищения направления р’ от накапливающихся 
погрешностей, дополнялась на каждой 5-й итерации условием р’ =У.Л"". Миними- 
зация методом сопряженных градиентов практически заканчивалась для 100 итераций. 
В сравнении с предыдущими методами он не существенно изменил невязку Ди’, что 
наглядно демонстрирует рис. 2, однако уменьшил максимальное расхождение най- 


денной и требуемой функции Т(г,2), которое составило пах|Т —Т.| =37К после тех 


же 100 итераций. 

Таким образом, применение классических экстремальных методов для минимиза- 
ции целевого функционала .Л в рассматриваемой задаче оптимизации с параболическим 
уравнением свидетельствуют о низкой эффективности данных методов. Как следует 
из рис. 1, 2, точное решение недостижимо за значительное число итераций. 
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Таблица 1 —- Минимизация целевого функционала „Л различными экстремаль- 
ными алгоритмами 


Ите- Метод Метод Метод Метод с регули- 
ра- наискорейшего монотонного сопряженных руемым направ- 
ЦИЯ 


спуска убывания градиентов лением спуска 


[9,8 _ 


Теперь рассмотрим решение задачи минимизации функционала /(и) прямым 
экстремальным методом с регулируемым направлением спуска [5], [6]: 
и (2) =и"(2)-Ь' 6" (2)У.Л(и";2) на 5,5" > 0, =0,1..... (12) 
Функция 0(2) регулирует сходимость и“ — и, согласно необходимому условию опти- 
мальности [5], [6] в форме 
У/Л(и‘;2) > 0 равномерно по 2, при и^ —>и.. (13) 
Алгоритм (12) с реализацией условия (13) позволяет решать бесконечномерные задачи 
оптимизации за конечное число итераций без преобразований управлений к конеч- 
номерным векторам. 


м, ЕДЖА ыь 


100 


Рисунок 2 — Сходимость функции и(2) к точному решению и. (2) 
при использовании метода (11) 


Проблемы практической реализации метода (12) заключаются в выборе функ- 
ций 0” (2) для удовлетворения условия (13). Для приблизительной реализации (13) 
использовалась идея шаблонных приближений на первой итерации [1], [5], [6]. Пусть 
первое приближение и'(2) является некоторой заранее известной функцией $ф(2) (шаб- 
лонная функция), градиент которой У./(ф;2) удовлетворяет условию (13) для А =1 
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(т.е. равномерно убывает после первой итерации). При этом из (12) можем найти 
5°(2)= ИО) ‚ УЛи?; 2) #09235. 
У.Л(и ;2) 
На следующих итерациях параметр &(2) не меняется. В данном методе от исследо- 
вателя потребуется сделать несколько пробных первых итераций для подбора подхо- 
дящей шаблонной функции Ф(2), удовлетворяющей (13). В рассматриваемой задаче 
было выбрано $ф(2) = 0,2и°(2). 

Здесь глубина спуска на каждой итерации вдоль выбранного направления ми- 
нимизации 6“ (2)У./(и';2) определяется числом Б* по методу (8). Если 9"(2)=1, то 
алгоритм (12) принимает вид алгоритма метода наискорейшего спуска (4). Функция 
6" (2) регулирует направление спуска и определяется на первой итерации из условия 
не более 20 % изменения функции и(2): 


0,2м° (2 
9" (2) = в 1 (14) 
У.Л(и ;2) 
Решение представлено в табл. 1, а на рис. 3 показана сходимость искомой функции 
к оптимальному значению и, для итераций 5, 10, 20. На последней итерации решение 


практически совпадает с точным тестовым значением и, при максимальном расхожде- 


= 0,3 К. 


нии найденной и оптимальной функций тах|Т —Т. 


Рисунок 3 — Сходимость функции и(2) к точному решению и. (2) 
при использовании метода (12) 


Имеется быстрая равномерная сходимость, в отличие от традиционных алгоритмов 
минимизации (4), (11). Из табл. 1 видно, что всего за 20 итераций целевой функцио- 
нал Л уменьшился в 12,2.10* раз, а невязка решений Ди в 190 раз, что демонстри- 
рует принципиально лучший результат, по сравнению с традиционными методами 
минимизации (4), (11). 

Данные результаты подтверждают высокую эффективность прямого экстре- 
мального метода (12) для решения задач. Метод (12) на несколько порядков лучше 
минимизирует целевой функционал по сравнению с традиционными методами при 
тех же вычислительных затратах. 
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Полученные результаты демонстрируют нецелесообразность, и даже некоррект- 
ность применения классических экстремальных методов для минимизации функционалов, 
связанных с процессами, описываемыми параболическими уравнениями. С другой 
стороны, экстремальные алгоритмы с регулируемым направлением спуска [5] типа 
метода (12) можно считать эффективными и рекомендовать их для решения подоб- 
ных задач. 
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Про зжн!сть екстремальних метод в в нескнченновим!рних задачах оптимзацй 

з параболчним рвнянням 

У стати розглядаеться зЫжнсть екстремальних методв у задач! оптимзацй системи з квазиинейним 
парабол!чним рвнянням. На приклад! нескинченновимрно! мимвацй квадратичного шльового функцщоналу 
показано, що традищйн! методи, що включають методи найшвидшого спуску, монотонного спадання 1 
зв’язаних градентв, збтаються не р1вномрно. Равномрно! з@1жност! дозволяе досягти прямий 
екстремальний метод з регульованим напрямом спуску. 


№А. Гоо4ат, У.К. То5букй 

АБош Сопуегоепсе о? Ехтете Мео4$ ТазК$ о Орйпи7айоп ш Зрасе Епез$пе$$ у РагаБоНс 
Еаца!7аНоп 

Сопуегоепсе оЁРехтете те фо4$ 15 ехаттей ш Ше ‘а5К оЁ орйпихайоп оЁ те зубет ул попПпеаг Бу 
рагабоЙс едиаПтайоп. [ 15 гойпе оп фе ехатр!е оЁ пшшитайоп ш расе епе5$пез$ оЁ диадгайс 
Бауше а зрес1а] ригрозе ЕапсНопа|, а та@@опа! те о4$, шсГадте Фе тефо4$ оЁ Фе теоа оЁ 1е 
гар! 1о\уегте, топоюпоиз 4есгеазе ап айеп4е4 ога 1еп{$, шеей по{ суету. АПо\з еуеп сопуегоепсе {0 
оБат а Ппе ехиете тео4 у Фе тапазе4 атесйоп о# ]о\уегпв. 
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